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1.  Contract  Description. 

The  focus  of  this  project  remains  to  investigate  a  wide 
range  of  phenomena  relating  to  energetic  materials  in  the  solid 
phase.  In  addition  it  is  proving  necessary  to  extend  the  range  of 
applicability  of  existing  techniques  and  to  develop  new 
methodology  to  enable  sufficiently  accurate  studies  at  an 
atomistic  level  on  these  little  studied  systems  (this  is  so  at 
least  in  the  level  of  traditional  solid  state  physics 
methodology).  The  need  to  develop  new  methods  also  necessitates  a 
few  simple  studies  on  prototype  systems  for  which  competing 
theoretical  calculations  are  available  and  also  for  which 
definitive  experimental  data  may  be  available.  All  theoretical 
studies  undertaken  or  proposed  are  to  be  of  the  ab  LDLtiS  type. 

In  the  case  of  solid  fuels  or  explosives  the  area  of  initial 
relevance  is  the  ini tiatiori’ of  the  combustion  or  detonation.  In 
the  real  world  this  may  be  by  means  of  a  strong  shock  wave 
impinging  on  the  system  and/or  by  means  of  elevated  temperatures. 
At  an  atomistic  level  this  means  one  looks  for  methods  to  lessen 
the  energy  required  to  produce  structural /chemical  modifications 
and  the  means  by  which  energy  is  transported  to  such  sites.  In 
practice  we  shall  see  that  such  modifications  are 
caused/accompanied  by  localization  of  energetically  excited 
states.  Such  localization  may  occur  either  at  defects  or 
impurities  or  may  spontaneously  occur  when  an  ’exciton’  causes 
sufficient  distortion  of  the  local  lattice  so  as  to  localize  in 
space  its  energy.  A  complication  of  this  study  is  that  the 
interesting  behavior  occurs  under  conditions  of  extreme 
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temperature  or  pressure.  Most  traditional  solid  state  or  quantum 
chemical  approaches  to  electronic  structure  determination  are  T=0 
type  theories  even  when  elevated  pressures  are  considered. 

After  combustion  is  initiated  it  is  essential  to  understand 
how  it  is  sustained  as  well  as  how  it  is  quenched.  The  process  of 
sustaining  combustion  may  well  be  autocatalytic  in  nature.  It 
seems  clear  that  the  portions  of  the  sample  near  the  point  of 
combustion  will  see  to  some  extent  at  least  the  same  extreme 
conditions  as  during  initiation.  The  modification  of  the  system 
by  such  extreme  conditions  therefore  takes  on  great  importance, 
and  our  early  efforts  have  been  directed  along  such  lines. 

The  ultimate  issue  of  stability  may  well  be  determined  by 
impurities  or  defects.  It  is  proposed  to  determine  what  the 
effect  of  defects  and  impurities  are  on  all  possible  properties 
of  these  molecular  systems  ^hat  may  relate  to 

combust ion /detonation.  It  is  clear  that  studies  on  both  surfaces 
and  bulk  systems  are  essential. 

2.  Short  List  of  Unsolved  Scientific  Questions  Related  to 
Molecular  Solids. 

If  one  excepts  the  solid  rare  gasses  and  to  a  lesser  extent 
solid  CH^,  and  such  essentially  one  dimensional  systems  as  the 
various  polymers,  few  if  any  atomistic  theoretical  studies  have 
been  successfully  performed  for  molecular  solids.  The  few 
available  studies  ignore,  other  than  for  several  solid  rare 
gasses,  systems  under  extremes  of  temperature  or  pressure.  At  the 
very  least,  the  studies  supported  by  this  contract  are  to  provide 
a  substantial  beginning  in  the  determination  of  such  fundamental 


properties  of  three  dimensional  molecular  solids  as  the  charge 
density,  the  electronic  density  o-f  states,  the  size  and  nature  of 
the  fundamental  electronic  excitations  and  possibly  atomic 
distortions  of  the  basic  unit  cell  upon  electronic  excitation  or 
under  exposure  to  elevated  pressure. 

In  order  to  fully  appreciate  the  difficulties  consider  the  case 
of  solid  CH^,  the  best  understood  solid  of  this  type.  Crystalline 

CH.  is  a  van  der  Waals  solid  in  which  the  carbon  atoms  form  in  a 
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close  packed  fee  lattice.  The  H’s  however  form  in  tetrahedral 
symmetry  about  each  C  atom  but  do  not  achieve  any  long  range 
order.  Thus  in  the  ground  state  the  basic  CH^  unit  is  free  to 
rotate  and  in  many  respects  behaves  as  a  spherical  molecule.  Does 
this  simple  behavior  extend  to  the  excited  configurations  as 
well? 

It  is  clear  that  littl'e  definitive  is  known  about  solid 
CH^,  and  it  must  be  conceded  that  even  less  is  known  about  a 
more  interesting  system,  solid  NO.  Nonetheless,  from  what  is 
known  it  is  quite  clear  that  there  are  definitely  peculiar 
properties  associated  with  the  geometry  and  electronic  structure 
of  this  solid.  The  crystal lographic  unit  is  <NO)^.  This  unit  does 
not  appear  to  have  a  strong  chemical  identity  however,  but  rather 
is  composed  of  two  N^O^  units.  The  evidence  for  this  is  the 
following:  the  separation  of  the  two  N^O^  units  is  large  and  the 
magnetic  properties  of  the  unit  <N0)^  are  similar  to  the  N^O^ 
unit,  being  diamagnetic.  This  is  in  some  ways  puzzling  in  that 
within  a  N202  unit  separation  of  the  two  NO  molecules  is  far 

greater  than  in  free  space,  making  strong  interactions  of  the  two 
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units  unlikely.  Furthermore,  the  individual  NO  molecule  is  known 
to  be  paramagnetic.  The  -fact  that  the  individual  N^O^  unit  is 
diamagnetic  indicates  some  real  interaction  between  the  two  NO 
units  must  exist  despite  the  large  separation.  This  all  raises 
important  questions  regarding  solid  NO.  What  is  the  nature  o-f  the 
fundamental  electronic  excitations  of  the  solid  form  and  how  do 
they  relate  to  those  in  the  free  NO  molecule?  Is  the  basic  unit 
able  to  maintain  its  geometry  under  pressure  and  if  not  then  what 
modifications  occur,  and  at  what  type  of  pressure? 

3.  Scientific  Approach  and  Recent  Progress  in  its  Implementation. 

Two  complementary  approaches  are  to  be  used:  the  first  of 
which  is  applicable  to  bulk  systems  and  to  ordered  surfaces.  This 
is  traditional  energy  band  theory.  In  its  current  implementation, 
we  employ  the  Linear  Combination  of  Molecular  Orbitals  method 
(LCMO)  in  a  self  consistent  Hartree-Fock  (SCFHF)  approach.  This 
level  of  theory  is  augmented  by  a  direct  application  of  (nearly) 
degenerate  perturbation  theory  including  corrections  to  the 
second  order,  to  enable  us  to  include  relaxation  corrections, 
correlation  corrections  and  the  formation  of  localized  electron- 
hole  pair  excited  states  (excitons).  Extensive  calculations  of 
this  type  have  been  performed  for  solid  CH^  and  the  results 
prepared  for  publication.  As  a  substantial  detailed  manuscript 
for  this  exists,  I  shall  not  repeat  the  lessons  learned  in  this 
report . 

The  second  approach  is  to  model  systems  by  finite  cluster 
techniques.  These  methods  are  of  particular  importance  for 
phenomena  in  which  substantial  departure  from  translational 


invariance  are  expected.  This  would  definitely  include  the  case 
of  lattice  distortion  accompaning  an  excitation  or  around  a 
defect  or  impurity.  In  the  implementation  of  the  finite  cluster 
method,  we  treat  all  electrons  and  nuclei  inside  the  cluster  in 
an  Unrestricted-Hartree-Fock  <UHF)  way  and  augment  the  UHF 
results  in  all  cases  with  a  perturbative  calculation  of  the 
correlation  properties.  In  the  case  of  near  degeneracies  we  also 
perform  a  nearly  degenerate  perturbation  treatment  as  well.  This 
is  of  value  in  that  by  such  means  we  eliminate  the  broken 
symmetries  of  the  UHF  method  and  recover  Multiplet  Structure.  In 
all  cases  we  solve  for  the  basic  UHF  state  in  a  self  consistent 
manner,  and  thus  spectral  detail  is  obtained  by  taking 
differences  of  total  energy, rather  than  from  some  more 
approximate  way. 

Substantial  effort  had  gone  into  determining  the  accuracy 
needed  in  perturbative  treatments  in  order  to  predict  spectra. 

Dr.  Beck  has  for  example  studied  first  the  Methane  molecule. 

Using  an  enlarged  basis  set  including  f— orbitals,  he  has 
succeeded  in  obtaining  90%  of  the  total  correlation  energy  in 
both  ground  and  excited  state.  In  doing  this  one  computes  to  high 
accuracy  and  great  expense  the  correlaton  energies  of  the  cores 
and  other  levels  not  involved  in  the  electronic  transition.  These 
correlation  contributions  largely  cancel  upon  electronic 
excitation  except  for  those  few  terms  associated  with  the 
ruptured  bonds.  We  have  studied  the  use  of  smaller  correlation 
bases,  which  primarily  describe  the  electrons  which  participate 
most  in  the  excitation.  Our  MBPT  code  was  organized  to  produce 


intermediate  output  which  -Facilitated  our  being  able  to  determine 

this  answer.  A  test  o-f  the  restricted  basis  was  made.  The  total 

correlation  energy  is  only  about  50%  o-f  the  experimental  value, 

but  the  change  in  total  energy  upon  excitation  is  essentially 

unchanged.  This  technique  is  extended  to  solid  CH^.  In  the  case 

o-f  a  bulk  cluster  we  use  13  CH.  units  -for  a  net  of  65  atoms.  All 
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electrons  are  treated  here.  The  excitation  energy  is  computed  and 

reexamined  upon  distortion  of  the  base  Methane  unit.  The 

fundamental  exciton  is  stable  against  ionization  by  an  amount 

greater  than  1  eV.  and  a  distortion  is  further  stabalized  by 

in  excess  of  1.5  eV.  This  is  consistent  with  the  earlier  energy 

band  theory  results  as  well. 

A  second  cluster  of  9  methane  molecules  was  treated  to 

simulate  a  surface  situation.  Again  the  simple  exciton  in  the 

ideal  T^  geometry  was  found, to  be  stable  against  ionization,  and 

had  an  excitation  energy  close  in  value  to  the  free  molecule. 

Distortions  in  response  to  this  excitation  may  be  in  either  a  C2v 

or  a  D_.  mode.  Both  are  considered.  The  lowest  excited  state 
Zn 

appears  to  have  a  C2v  geometry  in  which  the  elongated  bond 

lies  into  the  solid.  This  is  relevant  in  the  extreme  to  the 

present  study.  The  presence  of  an  electrical  excitation  at  the 

surface  of  solid  Methane  causes  an  exciton  which  spontaneously 

lowers  the  symmetry  of  the  system  and  partially  predi ssociates 

the  excited  CH.  molecule  into  a  CH_  unit  and  an  H„  unit.  Both 
4  2  2 

these  units  are  inherently  reactive.  A  report  of  these  studies 
should  soon  be  prepared  for  publication. 


5.  Publications  and  Related  Documentation. 
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